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Abstract 

The interplay between the electron repulsion U and the Jahn- Teller electron- 
phonon interation -Elr is studied with a large d model for the ferromag- 
netic state of the manganese oxides. These two interactions collaborate to 
induce the local isospin (orbital) moments and reduce the bandwidth B. Es- 
pecially the retardation effect of the Jahn- Teller phonon with the frequency 
Q is effective to reduce B, but the strong fi-dependence occurs even when the 
Coulombic interaction is dominating (U >> £'lr ) as long as £'lr > f^- The 
phonon spectrum consists of two components, i.e., the temperature indepen- 
dent sharp peak at a; = = ^1[{U + AE^^)/!!]^^^ and that corresponding to 
the Kondo peak. These results compared with the experiments suggest that 

^ < Elr < U in the metallic manganese oxides. 
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Since the discovery of colossal magnetoresistance (CMR) the manganese oxides 
Rei-xAj-MnOs ( Re: rare earth metal ion, A: divalent metal ion ) have been attracting 
intensive interest recently One of the controversial issues is the role of orbital degener- 
acy of Cg-electrons. In the conventional models of these compounds the strong Hund's 
coupling is considered to be of primary importance and the orbital degeneracy has been 
often neglected. Recently, however, Millis et al. P,^] correctly pointed out that the Hund's 
coupling alone is not enough to explain the CMR together with the insulating temperature 
dependence of the resistivity p(T) above the ferromagnetic transition temperature Tc. The 
additional coupling they proposed to be important is the Jahn- Teller electron-phonon cou- 
pling, which lifts the double degeneracy of Cg orbitals and gives rise to the polaronic effect. 
The small polaron formation above leads to the insulating p(T). In contrast with the in- 
sulating phase, the spin alignment below will enlarge the band width B as B oc cos{6ij/2) 
{Oij is the angle between the two neighboring spins Si and Sj), and the Jahn- Teller interac- 
tion enters the weak coupling regime. This scenario seems to be consistent with the large 
isotope effect on T^. by replacing O^^ by O^*^ |^. 

The strong correlation models of these compounds have been also studied by several 
authors especially for the undoped cases @JT^. In these models the Jahn- Teller interaction 
is considered to be smaller than the Coulomb interactions, and the effective Hamiltonian for 
the spin-orbital coupled system and its phase diagram have been clarified. It is rather natural 
to assume the strong electron correlation because the strong Hund's coupling originates from 
the strong electron-electron interactions ( ~ 5eV ). However, it is not trivial if the electron- 
electron interaction continues to be strong even in the effective Hamiltonian describing the 
low energy physics after the screening by oxygen orbitals and conduction electrons, and 
it still remains the controversial issue whether the electron-electron interation and/or the 
Jahn- Teller coupling are in the strong coupling regime or not. 

Experimentally there are several anomalous features which cannot be explained by the 
weak Jahn- Teller coupling described above even in the low temperature ferromagnetic phase 
below Tc. 
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[a] In the neutron scattering experiment no temperature dependent phonon modes have 
been observed [|1^. The recent Raman scattering experiment also shows that the Jahn- 
Teller phonons ( especially their frequencies ) are temperature independent and insensitive 
to the ferromagnetic transition at ||1 2|| . 

[b] The photo-emission spectra show a small discontinuity at the Fermi edge even at T << T^, 
which suggests some interactions still remain strong there UlSj . 

[c] The optical conductivity ct{uj) at T << Tc is composed of two components, i.e., the 
narrow Drude peak {u < 0.02eV) and the broad incoherent component extending up to 
uj ~ leV fl^. The Drude weight is very small, which seems to be consistent with the 
photo-emission spectra. 

[d] The low temperature resistivity p{T) can be fitted by 

piT)=po + AT', (1) 



where A is a large constant of the order of 500/if2cm/K^ |]T5[, again suggesting the strong 
electron correlation. 

[e] Contrary to the case of resistivity [d] , the coefficient of T-linear specific heat is very small 
with 7 = 2mJ/K, which violates the Kadowaki- Woods law for these compounds [|r^ . 



Although [e] is difficult to reconcile with [a]-[d], we consider the latter as the evidences 
for the strong coupling even at T <^ Tc. Because the spins are perfectly aligned at T ^ T^, 
the only remaining degrees of freedom are the orbital ones. In this paper, we study a large d 
model for the ferromagnetic state including both the electron-electron interaction U and 



the Jahn- Teller coupling g. This is the generalization of Ref.[7] in two respects, (i) including 
the electron-electron interation, (ii) including the quantum fluctuations. Especially the 
latter is essential to describe the low temperature Fermi liquid state, which is described by 
the Kondo peak in the large d limit []T^. The strong electron-electron interaction with the 
reasonable magnitude of the Jahn- Teller coupling explains both the large isotope effect and 
[a]. Moreover, the features of strong correlation [b]-[d] are at least consistent with the large 
U picture although [e] still requires further studies, which we have not undertaken in this 



paper. 

We start with the Hubbard-Holstein model for the ferromagnetic state. 

i,j,a,a' i 

- 9T.Q^i^^i - + \ + Mn'Q')^ (2) 

i i 

where the spin indices a =],[ correspond to the orbital degrees of freedom as t= d^^^yi and 
1= d^^2_r2, and the real spins are assumed to be perfectly aligned. We consider only one 
Jahn- Teller displacement mode Qi for each site, while there are two modes Q2i,Q3i in the 
real perovskite structure |jl8[ . However this does not change the qualitative features obtained 
below. The Jahn- Teller mode is assumed to be an Einstein phonon with a frequency Q. The 
transfer integral t^"' depends on a pair of orbitals (a, a') and the hopping direction. These 
dependences lead to the various low lying orbital configurations, and thus they suppress the 
orbital orderings in the ferromagnetic state. Actually there are no experimental evidences 
for the orbital orderings, e.g., the anisotropics of the lattice constants and/or transport 
properties in the ferromagnetic state. Then we assume that no orbital ordering occurs 
down to the zero temperature due to the quantum fluctuations, and the transfer integral 
is assumed to be diagonal for simplicity, i.e., t^"' = t^Saa'- This means that the ground 
state is a Fermi liquid with two degenerate bands at the Fermi energy. In order to describe 
this Fermi liquid state, we employ the large-d approach where the model Eq.(2) is mapped 
to the impurity Anderson model with the self-consistent condition This approach has 



been apphed to the manganese oxides to study the Hund's coupling by Furukawa [|] and 
to study the additional Jahn- Teller coupling by Millis et al. [0. The action for the effective 
impurity model is given by 



S= f^dr dT'ci{r)G,\T-T')c^{r') 
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2 



{drQ{T)f + n''Q{ 



T 



(3) 



f/n|(r)n|(r) - ^g(r)(n|(r) - n^^r)) 

where Gq^{t — r') is the dynamical Weiss field representing the infiuence from the sur 
rounding sites. The self-consistency condition is that the on-site Green's function G{iuJn) = 



{Go{iu;n) ^ — ^ should be equal to the Hilbert transform of the density of states 

D{e) as 

D{e) 



de 
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iUn + A* — Ti{i0Jn) — £ 

We take the Lorentzian density of states D{e) — t/{n{e'^ + 1'^)) because there is no need to 
solve the self-consistency equation in this case, and the Weiss field is given by 



Gq {iujn) = + + it sign 



(5) 



The only unknown quantity is the chemical potential /j,, which is determined by the electron 
number. The determination of chemical potential requires some numerical calculations, but 
our discussion below is not sensitive to the location of the chemical potential. Then we take 
II as the parameter of our model, and the problem is now completely reduced to that of a 
single impurity. Now we introduce a Stratonovich-Hubbard (SH) field ^ (r) to represent the 
Coulomb interaction U. Then the action is given by 

S = J2{iu;ri + n + itsignu;n)cl{u;n)caM 



+ / dr 

Jo 



^ii.rf + ^mQi.r)? + ^'Q{Tf] + C(r)Q(r) - ^r) + ^g(r))K(r) - n,(r)) 



where A = (i^ + and C(t) is the test field to measure the phonon correlation function. 
At this stage the electron is coupled with the linear combination of the SH field ^ and the 
phonon Q as rjir) = ^(r) + ^Q(t). Then the effective action can be derived as 



+ E 



2A2 



+ -JJ-{C{i^n)v{-i<-^n)+r]{i<-^n)C{-i<-^n)) 



(7) 



where Q = D,^JUeG./U. The effective interaction Ues. is the sum of the Coulomb repulsion 

U and the lattice relaxation energy i?LR = g'^/{2MVL^) as Ues. = U + 4£^lr- The phonon 
Green's function d{iujn) — {Q{ iuJn)Q{—iuJn)) is given by 



where x^ii^n) = {v{'^^n)v{~'^^n)) is the orbital susceptibihty of the electronic system. The 
first term is the usual phonon Green's function with the renormalized frequency Q, which 
is higher than the bare Q. This Q does not depend on the electron response function and 
hence temperature independent. When U ^ -Elr, ^ is much larger than Q and the first 
term of Eq.(8) is irrelevant in the phonon frequency region u; ~ ~ O.leV. In the opposite 
limit U ^ -Elr, the renormalization of the phonon frequency is small. The second term 
depends on the temperature (Kohn anomaly) and shows the characteristic lineshape with 
the broadening. From the above consideration together with the experimental fact that no 
temperture dependence of the phonon frequency has been observed, we conclude that the 
Coulomb interaction U is larger than the lattice relaxation energy -Elr- 



Now let us study the electronic system in detail. Following Hamann |T9[, we can obtain 
the solution of the integral equation for the electron Green's function in the presence of the 
time-dependent field //(t), and hence for the effective action. The result is 
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77(r) tan ^ rjij) — + ti^tY) 
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In 



(9) 



dr' T]{Tf-r]{T'f l + r]{T'f' 
where P denotes the principal value. The Eq.(9) differs from Eqs.(3.40)-(3.43) of Ref. fT^ 
by the frequency dependent coefficient of rj{iuJn)Tj{—iuJn) due to the retardation effect of 
the phonons. A similar model without the Coulomb repulsion U has been studied by Yu- 



Anderson |20|, which corresponds to the limit U ^ 0, ^ oo, Vt^U = finite in Eq.(9). 
First, we consider the saddle point solution for the rj field assuming the static configuration 
ri{T) = rjQ. Then the action for the static solution becomes 



So = P 



eS. 



r/otan rjQ - + 7]^) 

TT \ / 



(10) 



When f/cff./ (^A) < 1) there is a single minimum at the origin, and the double minima appear 



for f/cff./(7rA) > 1. We are interested in the limit of strong correlation f/cff./(7rA) >> 1, and 
in this limit the tunneling events ( instantons) between the two minima at r] = ±i]q = ±^2a' 
play an essential role in producing the Kondo peak at the Fermi level. The third term in 
Eq.(9) gives the non-local interaction between the instantons along the time axis, which 
corresponds to Jz term when it is mapped to the Kondo problem. The fugacity z of the 
instanton, which corresponds to J±, can be estimated following ref. [l^. Let tq be the width 



of the instanton, i.e., hopping time, and the action for a single instanton is given by 

^inst.(ro) = rT9n^^ 9 ^^T^ (^^n)'?mst. ("^^n) + 7T"0^eff. " InTQ. (11) 

The last term can be written in terms of the derivative /i(r) = (i?7inst.(T)/rfT, which is the 
localized function around the center of the instanton ( r = 0) with the width of the order 
of To and the integral of /i(r) over r is 2r]Q. We assume h{T) = (?7o/To)e^''^'/'^". Then the 
summation over Un can be easily done and we obtain 

4i?LR ~ 1 

5'inst.(^o) = 7^2:;r"/(^^o) + qTqUcs. - In To, (12) 
where /(a) = a(a + 2)/(4(l + Q;)^). Minimizing S'inst.(''o) with respect to Tq, we obtain tq as 



To = 6/f/eff. for f/eff. » 6fi and tq = 3/U,s. + ^{^/U^^y + QE^^/{9?U,^) for U,^, « QQ. 
Introducing the dimensionless variables x = -Elr/^ and y = U/Ei^^, the reduction factor R 
for the instanton fugacity z = rj^^e'^ is given by 



R = 2x^l + Ay~\ for y » l/msix{x,x^) (13) 
R= + ^x2(4 + y), for ?/ << l/max(x, x^). (14) 

The main conclusion here is that when x > 1, i.e., -Elr > ^, the reduction factor R is 
proportional to x even in the limit U ^ -Elr- This is because the overlap of the phonon 
wavefunction enters the tunneling matrix element even when the Coulomb interaction is 
dominating. 

In summary, the fugacity of the instantons and the long range interaction between them 
can be mapped to the anisotropic Kondo model with the Kondo couplings given by 



where A^(0) is the density of states for the Kondo problem and the reduction factor R is 
estimated to be Eqs.(13) and (14). Solving the scaling equations for J± and Jz, an estimate 
of the Kondo temperature is obtained for large R{» 1) as 

RUcS.' 



Tk ~ t exp 



A 



(16) 



This Kondo temperature gives the effective bandwidth B in the large-c? model, and the 



effective mass enhancement is estimated as m*/m = t/Tx ~ exp 



. This strong mass 



A 

enhancement manifests itself in the physical quantities as follows. The orbital susceptibility 
Xri{i^n uj + i5) has a peak at = with the peak height of the order of i]q/Tk and 
with the width of the order of Tk- The specific heat coefficient 7 should be proportional to 
5^^, while the coefficient A of in Eq.(l) scales as A oc B^^. The ferromagnetic transition 
temperature is expected to be proportional to 5, and hence we expect a large isotope effect 
when R > 1. In other words, the large isotope effect does not rule out the large Coulomb 
interaction U . The edge of the photoemission spectrum at the Fermi energy and the Drude 
weight are reduced by the factor m/m*. As mentioned above, many experiments except 
for those of specific heat seem to suggest that the manganese oxides in their ferromagnetic 
state belong to the strongly correlated region. Then it is natural that the large isotope 
effect is observed. The band width is rather difficult to pin down accurately, but a rough 
estimate from A in eq.(l) is -B ~ O.leV, which is comparable to Vt. This suggests that 
the suppression is not extremely large. Another experimental fact is that the temperature 



independent Jahn- Teller phonon is observed at the frequency near the original one |]T2 
which means that U is larger or at least of the order of -Elr- Then we conclude that the 
relative magnitude of the energy scales is < -Elr < U for the ferromagnetic phase in the 
manganese oxides. The small specific heat coefficient necessitates the physical mechanisms 
which have not been included in our model. One possibility is the inter-site process like the 
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orbital singlet formation, and another possibility is that the proximate quantum criticality 
is affecting the Fermi liquid behavior. 

In summary, we have studied a model of manganese oxides in the ferromagnetic state 
taking into account both the Coulomb repulsion U and the Jahn-Teller interaction E'lr. 
These two interactions collaborate to induce the local orbital moments. In the strong cou- 
pling case, i.e., Ues. = U + 4£'lr >> t {t : transfer integral), the overlap of the phonon 
wavefunctions (~ e~^^^^^) enters the tunneling amplitude between the two minima for the 
orbital moment. The phonon spectral function consists of two parts, i.e., the sharp peak 
at the renormalized frequency Q = il^JU^s./U and the broad peak with the width of the 
order of the band width. These results can reconcile the large isotope effect on Tc and the 
apparent temperature independent phonon spectrum assuming U > £Jlr > ^■ 
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